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Summary. Using the model of particle acceleration in shocks 
in stellar winds developed by Biermann (1993) and Biermann 
and Cassinelli (1993) we demonstrate that the nonthermal 
radio emission spectrum from the cataclysmic variable bi- 
nary star GK Per can be interpreted consistently. We pro- 
pose that the radio spectrum results from hadronic secondary 
electrons/positrons, which, in the optically thin case, give a 
characteristic radio spectrum with a spectral index of -1-1/3 
below a maximum and, for a strong shock in a wind slow 
compared to the shock speed, of -0.7 above the maximum. 
We test the proposed concept for consistency with the radio 
data for OB and Wolf Rayet stars as well as radiosupernovae. 
There is some evidence that this process is important for the 
radio emission from active galactic nuclei, both in compact 
and extended emission. We argue that such a radio spectrum 
is a unique characteristic for hadronic interactions. 

Key words: Acceleration of particles. Accretion disks. El- 
ementary particles. Close binary stars. Cosmic Rays, Radio 
continuum: stars 



Introduction 

1.1. The origin of cosmic rays 

The origin of cosmic rays as observed near earth is still not 
fully settled. Recent work proposes that this origin can be 
attributed to three major sites, a) stellar explosions into the 
interstellar medium (ISM-supernovae), b) stellar explosions 
into a former stellar wind (wind-supernovae) , and c) powerful 
radio galaxy hot spots. The predictions were made in Biermann 
(1993, paper CR I), and in Rachen and Biermann (1993, paper 
UHE CR I). Tests on a) stellar winds were made in Biermann 
and Cassinelli (1993, paper CR II), on b) explosions into a ho- 
mogeneous interstellar medium in Biermann and Strom (1993, 
paper CR III), on c) the chemical abundance and spectrum 
against air shower data in Stanev et al. (1993, paper CR IV), 
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and on d) the ratio of wind-supernovae to ISM-supernovae 
and the implications for the chemical abundances of Helium 
and hydrogen in cosmic rays in Biermann et al. (1995). The 
chemical abundances and the predicted spectrum for the ex- 
tragalactic component (paper UHE CR I) were checked against 
air shower data in Rachen et al. (1993, paper UHE CR II). This 
latter check provided an independent control on the spectrum 
and chemical composition of the galactic component beyond 
the knee, near 5 10® GeV. A definitive check on the chemical 
composition across the knee will probably be possible with the 
MACRO-EASTOP experiment (Battistoni and Navarra, priv. 
comm.), combining muon data with air shower data. Clearly, 
any concept of how particles, nuclei and electrons get injected 
and accelerated has to account for the nonthermal radio emis- 
sion in other sites which are believed to be basically similar. 

1.2. The nova GK Per 

The nova GK Per, which exploded in 1901 and has been ob- 
served ever since with great care, is such a case. The radio 
spectrum of its shell has been observed (Reynolds and Cheva- 
lier 1984) to be a power law of spectral index ~ —2/3, which 
is close to the limiting spectrum of synchrotron emission from 
a shocked region in a stellar wind in the case that the shock 
is much faster than the wind speed (exactly —2/3). We note 
that this radio spectrum is quite unusual for novae (see, e.g., 
Hjellming 1990), since in other cases where radio emission has 
been detected, it can readily be understood as free-free emis- 
sion. This is similar to the case of OB and WR stars (CR 
II), where free-free radio emission is also very common, and 
non-thermal radio emission is clearly less frequent (Abbott et 
al. 1986, Bieging et al. 1989). For OB and WR stars we have 
suggested that shocks in winds accelerate particles to sufficient 
energies to give the observed non-thermal radio emission. In 
Nath & Biermann (1994) we have been able to check the con- 
cept of shocks running through the winds of OB stars and then 
hitting a surrounding shell with the observed 7-ray line and 
continuum data. 

In this paper we try to understand the radio emission from 
the shell of the nova GK Per in order to check the concepts 
described in the series of papers on the origin of cosmic rays, 
and develop them further. 
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The model which we have developed for strong shocks in 
stellar winds (sec, MacFarlanc & Cassiriolli 1989, CR I and CR 
II) has been described in some detail in Biermann (1995). The 
model concentrates on spherical Shockwaves, with a magnetic 
field orientation mostly parallel to the shock front, as is the 
case for a spherical shock advancing into an asymptotic Parker 
type wind. Wc note that in this case a north-south asymmetry 
can develop due to the possibility that the sign of the particle 
drifts differs in the two hemispheres (towards the pole in one 
hemisphere and towards the equator in the other, for instance). 
We propose this effect to explain the asymmetry in the radio 
emission in the radio shell of GK Per. 

Therefore we explore in this paper the idea that the non- 
thermal radio emission from the shell produced by the nova 
GK Per is also due to a strong shock initiated in a wind by 
the outburst of 1901. For the wind we assume a Parker-type 
spiral magnetic field (see, e.g., Jokipii et al. 1977). 

2. 

The low frequency turnover in the radio spectrum 

2.1. The data for GK Per 

The spectrum of GK Per has been published by Seaquist et al. 
(1989). We have reanalyzed the WSRT data at 49 and 92 cm, 
making full use of the redundancy technique to achieve high 
dynamic range and thereby reduce the interfering effects of 
background sources. This results in an upward revision of the 
flux densities Sv at both wavelengths, and smaller errors. The 
flux density scale is based upon that of Baars et al. (1977), 
using 3C48 as a primary, and 3C295 as a secondary calibrator. 
The flux densities are 8.7 ± 0.5, 20.6 ± 1.6, 29 ±3, 33 ± 10, 
38 ± 6, and < 40 mjy at the frequencies 4.86, 1.49, 0.608, 
0.408, 0.327, and 0.151 GHz, respectively. The improvement 
over previous data is the Westerbork points at 0.608 and 0.327 
GHz. Our spectrum is shown in Fig. 1. It is clear that the high 
frequency spectrum turns over, gradually, somewhere below 1 
GHz, suggesting either absorption or a low energy cutoff in 
the electron spectrum. 

When we fit a powerlaw through the three high frequency 
data points, we obtain 



= (13.4 ±1.3 mJy)(^^)-0-^3±0.08^ (1) 

where we have used the geometric mean between the high 
frequencies as reference. The deviation of the observed flux 
densities from this fitted powerlaw is 2 to 3 a for each indi- 
vidual frequency only, but at all four lower frequencies this 
deviation is in the same direction, malcing it of interest to 
pursue; to have such a deviation occurring systematically is 
signiflcant. Since systematic errors enter the individual errors, 
it is difficult to measure the significance of this systematic 
difference, but assuming the systematic errors to be fully ac- 
counted for in each individual error determination, then the 
combined significance is high, much higher than 3 a. In the 
following this observation is the key to the arguments we will 
be presenting. 



Fig. 1. The newly derived radio spectrum of GK Per using WSRT 
data in the redundancy mode. These newly derived data points are 
emphasized in this figure which is otherwise analoguous to the spec- 
trum shown in Seaquist et al. (1989) 
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2.2. Interpretation 

What is the physical origin for this downturn at low frequency? 
It is very difficult and quite unconvincing to attribute the 
cutoff to free-free absorption, or synchrotron self- absorption 
(Seaquist et al. 1989). We note that Seaquist et al. (1989) 
already mentioned the possibility of a low energy cutoff in the 
energetic electron distribution. A low energy cutoff in the elec- 
tron distribution is a natural consequence of miclcus-uuclcus 
interaction betwen cosmic ray particles and interstellar mat- 
ter, primarily proton-proton (pp) interaction; this is due to the 
large rest mass of charged pions which are created in such col- 
lisons (neutral pions decay into photons, of course), and which 
then decay further, ultimately leaving neutrinos of various fla- 
vors and electrons and positrons. This topic of the injection of 
secondaries from energetic particle interaction with matter has 
been of classical interest (Hayakawa k. Okuda 1962; Ginzburg 
& Syrovatskii 1964; Ekspong et al. 1966; Stecker 1971; and 
many others; recent references are Klcmcns 1987, Bicrmann & 
Strittmatter 1987, Sikora et al. 1987 and Mannheim 1993a, b). 

The kinematics of the creation of a new particle like one 
or several charged pions in the center of mass frame of the 
encounter of the two protons reduce the energy available for 
the pion when the excess energy comes down to the scale of 
the proton rest mass, and so we expect the electron/positron 
distribution to deviate slightly from the ultimate power law 
below energies of a few GeV. Such a change is indeed observed 
in proper Monte-Carlo calculations (Daum & Biermann, in 
prep.). The peak of the production spectrum, when multi- 
plied by energy to find the relevant maximum for synchrotron 
emission, is near 100 MeV: 7e,mm 

« 200. Thus the final elec- 
tron/positron spectrum is the well known power law above 
a particle energy of order 5 GeV, but shows a low energy 
turnover and cutoff at about 100 McV. Between 100 McV 
and about 5 GeV the spectrum is somewhat flatter than the 
ultimate high energy spectrum. 

Hence, the expected radio spectrum is a power law Sv ~ 
v^^^^ in frequency v below 

Vmax = 7.2 10^ Hz O^^f sin 6'cone Sq.s, (2) 

where we have referred the magnetic fleld at the actual location 
of the shock to the magnetic field at the equator at the ref- 
erence radius by B{shock) = 1.1 10~ 

sin ^corie -^0.5 from the 
assumption of a Parker wind topology; Bq.S is the strength 
of the magnetic field at the equator at the reference radius of 
10^^ cm in units of 3 Gauss. 9cone is the cone opening angle 
of the emission region. 

The low frequency turnover is in the range 0.3 — 1.4 GHz, 
and so this suggests a magnetic field strength of about 1.3 to 
5.8 Gauss at the reference radius, and more if the cone angle of 
the observed emission is small. Here the latitude dependence 
of the magnetic fleld strength enters only linearly, and so the 
magnetic field implied here is the lower limit for the equatorial 
magnetic field strength. We thus estimate from the upper as 
well as lower limit the magnetic field to be of order 1 Gauss at 
the reference radius. 

Above the frequency for the maximum in emission (corre- 
sponding to 7e,mi?i) the radio spectrum has the approximate 
form Si/ ~ f "'^j for Ui/V\y ^ 1, but somewhat steeper than 
this for smaller shock velocities, with some steepening for the 



frequency range a factor of order 3000 above the frequency of 
maximum emission. Here Ui is the shock velocity in the wind 
frame, and W is the wind velocity. We emphasize that a) this 
spectrum is optically thin, and b) the spectrum is very broad 
at the maximum of the flux. 

We have used here the Fermi model for the multiplicity of 
the secondary pions; there is a weak modification for the spec- 
trum of the secondaries when one considers the full kinematics 
of the encounter of the pp-coUison, measured multiplicities and 
the thresholds (Daum & Biermann, in prep.). 

This spectrum is valid for primary proton spectra of E~'^^^ 
which have been derived for strong shocks which arc much 
faster than the wind; for steeper proton spectra such as ex- 
pected for shocks which are slow compared to the wind the 
resulting secondary spectrum is also steeper, and so we cam 
reconcile the observed radio spectrum and its error range with 
this suggestion. 

This spectral behaviour of two well-specified power laws 
around a peak is thus a characteristic signature of hadronic 
interactions, and provides a clear statement: Acceleration of 
nuclei, including protons, is clearly indicated by such a spec- 
trum. Then there ought to be corresponding 7-ray emission, as 
well as neutrino emission, which, however, usually is too weak 
to be detected. We note that the threshold for new particle pro- 
duction is a common signature for hadronic interactions, see, 
e.g., also Mannheim (1993a) leading quite generally to breaks 
between two different power laws. Whether such processes also 
help understand 7-ray bursts (e.g.. Band et al. (1993)), which 
also seem to show such nonthermal spectra, is open at present. 



3. 

Injection of secondary particles 
3.1. Basic balance of injection and loss 

In this main section of the paper we will derive the equivalent 
injection efficiency of secondary particles from p-p interactions 
in order to see whether we can explain the radio data and their 
implications for i) OB and WR stars, ii) radio supernovae, 
and most importantly, iii) GK Per. The basic question is to 
ask is whether we can explain the rather different injection 
efficiencies as due to primary injection from the thermal pool 
in some cases, and as secondary injection in all other cases; 
the test is the quantitative estimate of the efficiency derived 
and, also, obviously the spectrum of the particles. 

We use the results and notation from CR II and define 

Vppul = ^ Cpt-^^^Ekindt. (3) 

with t := pp/{mpc); herepp is the momentum of a proton, E^.^„ 
its kinetic energy, Ui the shock speed, and rjp an efficiency 
of converting shock kinetic energy to an energetic particle 
population. The integration limits ta and tf, are from the 
beginning of the energetic particle population close to the 
thermal MaxwoUian tail, to the maximum achievable; as long 
as <C 1 and if, 3> 1 the integral is independent of their 
exact values. Below we will refer to the equivalent efficiency 
for electrons rje directly obtained from observations, and our 
model for the emission region (see paper CR II). The density 
p above is the density in the wind of speed Vw and so we have 
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Cp = 26.9 7?p,_i( 



M_ 



Vw,- 



) ul. 



(4) 



'14 



The units for Cp are cm '^,10 ^ M, 



^yr for the mass loss, 10^^ 



cm for the radius, 0.1 for rjp, and 10~ c for the velocities. Then 
interactions between energetic protons and thermal protons in 
the wind produce secondary electrons and positrons as a rate 
per unit volume 



J 



7p "^^^ {apf{-fp))cnd^p, 



(5) 



where f{^tp) is the multiplicity. The combination of cross sec- 
tion and multiplicity (crp /(7p)) we approximate as a constant 
of value 40 10"^"^ /a cm ; there are various calculations of 
this in the classical literature, e.g. Gould & Burbidge (1968), 
Stecker (1971), Stephens & Badhwar (1981); we take the par- 
ticular value from a calculation by Niomcycr (1995), which 
improves upon Gould & Burbidge (1968) and agrees with 
Stecker (1971), and obtain an estimate for of 1.4. As a 
lower threshold we use the approximation 7p = 1.4 from the 
threshold condition for the Deltarresonance (derived from the 
energy where the cross section falls to half its peak; value) . 

The rate of generation of secondary electrons/positrons is 
then 
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Here we have made no distinction between hydrogen and other 
nuclei in the gas of the stellar wind, but have just counted all 
nucleons. 

We consider the shock frame as it advances. Thou wc also 
have losses, first of all dilution losses - from the continuity 
equation in spherical coordinates - with the rate 



further acceleration can not change this number, but modifies 
the spectrum. This can thus be written as 



9 ^ -13/9 
^ - -^^eP^epl ■ 



(10) 



Replacing the radius r with the corresponding time depen- 
dence we have 



ri4 = 610 ^(f/i,_2 + Viy,_2)f, 
we have with 

A =1.1810"r?p,_i/zi 

^ M-5 ^2 1 1 



(11) 



(12) 



the balance equation 



dF _ 
dt ~ 

where 



A 



e = 2 + 



c 



1 + Vw/Ui ■ 



(13) 
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Obviously, we have 2 < e < 3; the limit of Ui/Vw 2> 1 
and C = corresponds to e = 3. The integration yields 



^=3^i«T)^-(T)^)' 

for e ^ 3 and 



(15) 



Vw + Ui 



(7) 



and the losses downstream from the diffusive Fermi process 
(see Drury 1983) of 



(8) 



in the limit, used hero throughout, that the shock velocity 
ratio U1/U2 = 4; this is implied by our assumption that the 
shock is highly supersonic. The term C, describes the resupply 
of secondary particles from other sources further upstream; 
C = 1 denotes full resupply, i.e., no net loss, and C = is for 
no resupply. 

These processes apply to the number of energetic secon- 
daries, which is given by 



Ce 



-22/9 , 
7ep 07ep, 



(9) 



where the integration limits are 7epi ~ 200 and 7ep2 ^ 7epl- 
We concentrate on the number of secondary particles, since 



F = 4 In-, 

to 



(16) 



for e = 3. 

Now translating this expression back into a radial depen- 
dence we obtain 



= 2.55 10" 



' 14,0 
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and a corresponding expression in the two cases of e 7^ 3 and 
= 3. 

From this expression we can derive an equivalent value 
of r]*p which represents the energy of this population which 

we derive from observations ignoring the lower energy cutoff. 
The true energy content is lower by about a factor of 15 
because of the low energy cutoff (for lep,eff = 200). Using 
lep,eff = 7ep,l = 200 and = 1.4 this corresponds to an 
equivalent efficiency for electron/positron injection of 



5 



ntp =3.510' 
Vp-1 ( 



ri4,o r 
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r ' 3 — e 



(18) 
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Ui,-2 1 + vvM' 



and again a corresponding expression in the two Ccises for the 
value of e. This is, of course, only a good approximation as 
long as 7ep,e// » 1- 

The newly injected electrons and positrons are also subject 
to local energy gain and losses, which flatten their spectra to 
an index of —7/3 with the same low energy cutoff (sec Drury 
1983, section 2.3.2, with an added integration over the various 
source energies). As a consequence the numerical factors above 
change only slightly. 

In the case that particle interaction is strong initially, say 
when the shock is close to the material of the accretion disk, 
then the forrrmlae above change in an obvious way. In the 
balance equation the term corresponding to the newly created 
particles is zero and only an initial population exists which is 
slowly being diminished. There is an initial equivalent efficiency 
and thereafter it changes radially as 



, 7ep,e//(T-) 13/9 To 

°^ 200 ' ^ - 



(19) 



In the case that there is a continuous or quasi-periodic 
resupply of energy to the population of secondary particles 
^ ~ 1, this moans simply that ri*p = ri*p „ where »7*p.o can 
approach the energy density of primary particles riep,o ~ rjp. 
A limit is ^ep,eff{''') = const « 7ep,i. 

The induced numerical value for rjep has to be compared 
with data. Therefore we have to ask three questions: First, 
can we explain the estimated value of rje of order 10~^ for 
low velocity shocks in OB and WR stars, can we interpret the 
value of about 0.1 in SN-wind shocks, and can we at the same 
time interpret the value required by the observations of GK 
Per, of order nearly unity? This will then also give information 
on what to predict in secondary /primary ratio for high energy 
electrons. In order to do this we need to consider first the 
injection of protons. 



3.2. Injection of energetic protons 
3.2.1. The perpendicular case 

Therefore we have to discuss the injection for protons (and nu- 
clei), since we need to derive a physical theory for the value of 
To, where the protons start getting injected cfRcicntly into the 
acceleration process (see Bell 1978a, b, Drury 1983). Consider 
a proton in the post-shock thermal tail of the Maxwellian 
distribution. In many papers Ellison and collaborators have 
shown that injection from the thermal tail downstream is pos- 
sible (see, e.g., Ellison et al. 1990, Jones & Ellison 1991). We 
assume that a substantial population of such particles is built 
up. 

We first consider the case where the unperturbed magnetic 
field is nearly parallel to the shock front; this case is valid over 
most of 47r for a spherical shock going through a Parker wind 
in its limiting configuration. In a Parker wind there is a polar 
cone region, where the magnetic field configuration is radial, 
and thus is perpendicular to the shock front. We consider this 
second case in the following subsection. 



Now we ask the question: Can the injection of particles 
be maintained to keep on producing the energetic particle 
population? Consider for the sake of an example the case of 
a proton: Such a proton is subject most of all to two-body 
encounters with the thermal electrons, also simply referred 
to as ionization losses. The slow and steady energy gain from 
Fermi acceleration including drifts is hindered when the energy 
loss during one acceleration cycle is equal to or larger than the 
energy gain per cycle. The rate of energy loss is given in, e.g., 
Nath &: Biermann (1993), and this rate has to be multiplied 
by the residence time in order to obtain the loss per cycle. 
In the subrolativistic regime the change in velocity /3c can be 
written as 



1.02 10 
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where we have used a factor of 1.3621 rup between mass density 
and hydrogen density and a factor of 1.181 between hydrogen 
density and electron density, valid for a fully ionized gas at 
cosmic abundances (Schmutzler 1987). 

The gain per cycle in Fermi acceleration can be written 
in the subrelativistic limit as {Af3) 



gain 



10~^C/i _2. Putting 
now gain and loss equal gives a lower limit for the radius, 
where shock acceleration can be sustained, Vo'- 



r„,u = 0.10(i^^) 1 
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Putting this back into the expression for the equivalent effi- 
ciency rjep we obtain for e ^ 3 



V*ep =4.810-^(^)^-2 
r 
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and correspondingly for e = 3, where we have also included 
now the small correction for secondary acceleration of the 
secondary particles. 

3.2.2. The parallel case 

A contrast to the picture above is injection in the polar cone 
region, because there the magnetic field is predominantly radial 
already in its unperturbed stage. In the injection argument this 
means that we have to replace the residence time obtained from 
the concept of fast turbulent convection with a residence time 
derived from a turbulent wave field of waves parallel to the 
magnetic field. 

The residence time (see, e.g., Drury 1983) is given by 



Tres 



4k 
Ui0c' 



(23) 



where k is the scattering coefficient, /3 is the particle velocity in 
units of the speed of light c, and Ui is again the shock velocity. 
In the limit of strong turbulence at all relevant wavenumbers 
the scattering coefficient is given by 
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(24) riep 



where B{r) is the radial magnetic field, dominant in the polar 
cone region region. We use here the nonrelativistic limit. Prom 
Parker's wind theory we have for the polar cone region, that 



B{r) = 310-^GaussBo.5-^ 

J /its 
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(25) 



where O is the angular rotation rate of the star, Vw is again 
the wind velocity, and Rs,ii is the radius of the star in units of 
10^^ cm. The magnetic field topology changes from a predomi- 
nantly radial configuration near the equator to an increasingly 
tangential configuration at the radius given by V\y/^- As a 
consequence the residence time can be written in this limit as 
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4.6 sec 
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(26) 



This time is so short because we use the extreme limit of 
the Larmor radius of the particles for the characteristic length 
of the transport. 

The loss against interaction with the post-shock hot elec- 
trons is the most effective interaction, and peaks at that veloc- 
ity for the particles when the protons have the same velocity 
as the heated thermal electrons. As before we use post-shock 
parameters for the loss arguments. Thus the interaction is max- 
imized at /3 = 0.23f7i^_2. The interaction leads then to a loss 
of momentum Af3 equal to the gain from the Lorontz trans- 
formation of a particle going through the shock and getting 
scattered back at the condition 
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which is independent of radius; this is because the quadratic 
dependence of density and magnetic field in the polar cone 
region cancel here. If this condition is met, then the injection 
of energetic protons is hindered, and secondaries cannot be 
produced. 



= 5.5 10" 
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suggesting that secondaries do not play a dominant role in 
the radio emission of OB and WR stars; the number is too 
small because it refers only to the polar cone region. Inverse 
Compton losses further diminish any role of these secondaries. 

The critical Alfvenic Mach number for shocks to inject 
primary electrons into an acceleration process in a shock is in 
the range of 3 to 30 - the value is highly uncertain at present 
(see paper CR III). We estimate the equatorial Alfvenic Mach 
number of shocks in the winds of OB and WR stars to be 
of order 0.3. Hence the angular range of electron injection is 
limited by 9cone ~ 0.01 to 0.1. Using then the implied conal 
area which participates in the radio emission, and the latitude 
dependence of the magnetic field of the Parker model wind, 
we arrive at the overall efficiency, i.e. that efficiency which we 
infer from observations ignoring the angular limitation. This 
number is then in the range 10~^ to 10~® as required by the 
observations. 

Our conclusion hero then is that primary injection in a polar 
cone can explain the low apparent efficiency of the nonthermal 
radio emission in OB and WR stars. We emphasize that we 
consider single stars here. Given sufficient flux density, this 
conal restriction may be observable with VLBI. 

3.3.2. Radio supernovae 

Second, we can ask what is expected for radio supernovae: 
Here we likely have e = 3, f7i,_2 = 3. The observed radio 
supernovae all except 1987A seem to have exploded into red 
giant winds with low velocity with ~ 0.01. It is easy to 

see that the radio emission from radio supernovae is fully con- 
sistent with all energetic electrons being primarily accelerated 
electrons which dominate by about two orders of magnitude 
or even more over the secondary electrons/positrons. However, 
when the supernova shock hits the surrounding shell that the 
steady wind often builds up over the lifetime of the star, then 
substantial secondary particle production takes place, which 
dominates over the earlier history of the shock running through 
the wind (see, e.g., Nath & Biermann 1994). We will discuss 
cosmic ray electrons elsewhere (paper CR VI, in prep.). 



3.3. Application to stars 
3.3.1. OB and WR stars 

Pirst, let us consider the normal OB and WR stars, where we 
have argued (CR II, and Biermann 1995) that rje is of order 

io-6±i. 

It is clear that for ordinary parameters for OB and WR 
stars the condition derived above, eq. (27), is never met. 
One can show analoguously that ionization losses and pp- 
interaction cannot readily limit particle acceleration in the 
polar cone region. 

This entails that injection of secondaries from pp-coUisions 
may start quite close to the star, where densities are high. 
The radial behaviour of r?ep is estimated from our expressions 
earlier with e w 0.5: 



3.3.3. GK Per 

Last, let us consider the case of the nova GK Per: Optical data 
on the emitting mass, usually identified with the ejected mass, 
(Pottasch 1959) give a mass of 7.0 10~^ Mq. We identify this 
with the mass in the shocked shell of the wind material, under 
the premise that the ejecta caused a shock to run through a 
pre-existing wind. With this notion we obtain for M_5/F^ _2 
a value of 0.22 /p. However, nova shells are not spherically sym- 
metric and so the spherically averaged mass in the wind shell 
is estimated to be approximately three times higher (Seaquist 
et al. 1989). We will use this higher value in the following, and 
use fp to account for any further differences. Using such an 
approach the uncertainty in To becomes so large, that we keep 
this radius in the expression and so to a fair approximation we 
have for GK Per 
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77ep =2.310" 



3-e e-2 
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/ Tep,e// .13/9 . ^ 

I 200 ' [/i,-3 • 



(29) 



The lowest reasonable limit for the injection radius is clearly 
the size of the binary system, or the inner region of the accre- 
tion disk, giving ro,i4 > 410~^ and ro,i4 > 410~^. With the 
full limitation of the observed radio emission to a cone and no 
allowance for clumping in the magnetic field we find that the 
efficiency can indeed be as high as required by the observa- 
tions, under the condition that the shock speed is smaller than 
the Alfven velocity at the equator. 

If we try to find solutions where there is appreciable sec- 
ondary production in the wind, we fail; the main injection, in 
this picture, has to arise from the inital environment, at the 
inner edge of the accretion disk. 

The source parameters derived here suggest that secon- 
daries from hadronic interactions, i.e. p-p interactions here, 
give an equivalent injection efficiency between ~ 10~^ and 1. 
All this implies, that 'yep,eff is approximately constant, which 
leads to a solution close to the C0.S6 

« const and maximal, 
and e very close to a value of 2. This in turn implies either a 
slow shock {Vw ^ or substantial resupply of secondaries 
from upstream (C ~ !)• 

Hence we have found a consistent, albeit speculative picture 
that explains the radio emission of GK Per as produced by 
secondary electrons/positrons, with a theory which also gives 
acceptable estimates for OB and WR stars, as well as radio 
supernovae, where the standard acceleration from the thermal 
pool is probably dominant, but restricted to the polar cone as 
an injector for OB and WR stars. 



Consequences for other sources and other 
wavelengths 

4.-1- Other stellar binary systems 

The model implies copious pp-interaction at the beginning of 
the shock travelling through the wind. A consequence is both a 
power law continuum from 7r° decay as well as an annihilation 
line from secondaries. In nova Muscac Chen et al. (1993) have 
observed a 511 kcV aimiliilation line, which they argue arises 
from a region near the inner edge of an accretion disk. In 
the same nova Goldwurm et al. (1992) describe observations 
of a 7-continuum with slope —2.38 ± 0.05, consistent with a 
primary spectrum of protons of slope near —7/3. 



4-2. Other radio sources with related characteristics 

There are a number of other radio sources which appear to 
show properties consistent with the picture which we propose. 
Secondaries have been discussed to explain other features in 
the spectra of radio sources (see, e.g., Klemens 1987, Biermann 
& Strittmatter 1987, Mannheim 1993a, 1993b). 

Klemens (1987) described the X-ray emission of the M87 jet 
successfully with the concept that secondaries from hadronic 
interaction give rise to synchrotron X-ray emission. 

The Galactic center source shows a spectrum which rises 
approximately as v^^^ through the radio range; it is not clear 



at this stage whether this spectral behaviour is due to optical 
thickness effects or due to another process such as discussed 
here. The Galactic center source has been discussed elsewhere 
(Falcke 1995), following our earlier modelling (Falcke et al. 
1993, Falcke & Biermann 1995, Falcke et al. 1995). 

The most prominent example among radio galaxies is Cyg 
A (Muxlow et al. 1988). The radio emission of the lobes shows 
a low frequency turnover near 200 MHz which then can be 
connected with an estimate of the relevant magnetic field 
to a minimum energy of the electron distribution of about 
71 — 350 ± 50. There is evidence that radio galaxies in general 
tend to flatten towards low frequency (Laing & Peacock 1980, 
Mangalam & Gopal-Krishna 1994), which may be related to 
a low energy cutoff in the electron distribution; the effect 
should be strong only at very low radio frequencies due to the 
estimated low magnetic field strengths in radio galaxies. 



Is there an alternative? 

Where is the main weakness of the model proposed? The key 
is clearly the argument that hadronic secondaries are able 
to account for the spectrum. Therefore, in this section we 
propose to ask whether primary injection and acceleration can 
also produce the pattern of radio polaxization, radio emission 
strength, radio structure and radio spectrum. 

The radio polarization has been argued to arise from the 
instability of a cosmic ray mediated shock (see paper CR I, 
CR II and Biermann 1995). The argument on secondary or 
primary injection is indifferent to this, since we use primary 
protons to drive the instability. 

The radio emission strength could be explained by invoking 
little or no magnetic clumping and just using the decrease of 
total flux from the geometry of the underlying magnetic field. 
Since the magnetic field strength in a Parker wind decreases 
as sin^ towards the poles, the total radio emission is reduced 
by a factor of one to two orders of magnitude, depending 
upon which cutoff angle we use. This entails a correspond- 
ing increase in the implied efficiency to about the maximum 
possible. The cutoff angle would still arise from the limit- 
ing condition on injection, and hence correspond to a lower 
magnetic field strength than used above. This runs into the 
problem that the efficiency may then have to be increased 
beyond any reasonable level. 

The radio structure would derive from the same argument 
as before, the drift asymmetry in the accelerated particle pop- 
ulation in the two hemispheres (see paper CR I, CR II and 
Biermann 1995). 

The radio spectrum would also be easy to accomodate in 
the concept that the explosion advances into the interstellar 
medium (paper CR HI); however, the spatial asymmetry is 
then difficult to interpret. 

Finally, for the low frequency cutoff we would require a new 
and rather difficult explanation, like foreground absorption, 
an extremely clumpy absorbing medium, or something else. 
Seaquist et al. (1989) discussed all likely possibilities, and 
found no other reasonable explanation. Therefore the radio 
observation of the low frequency turnover is the key observation 
on which everything rests. A very low frequency spectrum of 
GK Per would be decisive to settle the interesting question on 
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the role of hadronic secondaries, and thus on the presence of 
accelerated protons in a specific radio source. 



6. 

Conclusion 

We have tested a model proposed earlier to account for the 
spectrum, particle energies and chemical composition of cosmic 
rays with the observations of the nonthermal radio emission 

from the nova GK Per. Wc obtain a consistent picture implying 
a pre-existing strong stellar wind. Such a strong wind can 
strongly infiuence the evolution of a low mass binary like GK 
Per due to mass and angular momentum loss from the system 
(see, e.g., Kim et al. 1992, Iben 1995). 

We interpret the radio spectrum of GK Per as resulting 
from hadronic secondary electrons/positrons, which, in the 
optically thin case, give a characteristic radio spectrum with 
a spectral index of +1/3 below a maximum and, for a strong 
shock in a wind slow compared to the shock speed, of « —0.7 
above the maximum. Wo argue that such a radio spectrum is a 
unique characteristic for hadronic interactions. We can account 
for the observed radio emission, its spectrum including the low 
frequency cutofi:, its radio polarization properties, geometric 
shape and its luminosity. 

We have applied the concepts developed to OB and Wolf 
Rayet stars, as well as radiosupernovae, and have been able to 
verify that the efficiencies of injection implied by the observa- 
tions can be readily interpreted quantitatively, however, with 
primary injection. One consequence is that the radio emission 
in OB and Wolf Rayet stars is from a polar cone, just as in 
GK Per, but restricted much more in angle. 

We have shown, that the radio emission properties of the 
old nova GK Per can consistently be interpreted in the context 
of our proposed model for the origin of cosmic rays, and have 
made a small step in the direction of accounting for secon- 
daries from nucleus - nucleus interaction in the populations of 
energetic particles. Wo point out the key radio observation to 
confirm the model proposed, an accurate very low frequency 
spectrum of GK Per. 
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